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Abstract 

The concept of an optimal electron-phonon interaction spectral density as an 
Einstein spectrum which allows to describe all physical properties of a super- 
conductor in an optimal way is developed from Carbotte's original definition 
of an optimum spectrum. It is shown, using the borocarbides YM2B2C and 
LUN12B2C as examples, that such a concept is meaningful even for anisotropic 
systems. An Einstein spectrum is sufficient for clean-limit systems, a 2<5-peak 
spectrum is better suited for anisotropic systems with impurities. 



I. INTRODUCTION 

Conventional superconductors are well described by Eliashberg theoryi which treats su- 
perconductivity as a boson-exchange phenomenon. The dominant feature of this theory is 
the electron-phonon interaction spectral function a 2 F(u>) which can be determined from tun- 
neling experiments! or theoretically from band structure calculations. Using such an a 2 F{uo) 
within Eliashberg theory allows to reproduce the superconducting properties of a conven- 
tional superconductor within experimental accuracy and this established the phonons as the 
exchange boson between the two charge carriers building the Cooper pair in conventional 
superconductors. 

Concentrating on isotropic systems Carbottei developed the concept of an optimum 
spectrum based on earlier work of Leavens@ and Mitrovic and Carbotte.i Such a spectrum 
can be developed from a theorem which states that for a given strength A = J °° duj a 2 F(u) 
of the spectral density a 2 F{uj) the best shape that will maximize the critical temperature 
T c is a delta function spectrum 

a?F(u) opb = AS[u-u*(jj*)], (1) 

with the delta function placed at the frequency u*(fi*) at which the functional derivative 
5T c /5a 2 F(u>) displays its maximum for a fixed value of the Coulomb pseudopotential fi*. 
Carbottei extended this concept to encompass other physical properties such as 2A(0) / A;_bT c , 
the zero temperature gap A(0) to T c ratio, and a number of others. This concept establishes 
that a relation 

X = Ax{ii*) (2) 

always exists, where X stands for T c , 2A(0)/fceT c , etc. and is a universal number 

determined from Eliashberg theory for each property X and which varies only slightly with 

In essence the optimum spectrum gives information about the phonon frequency impor- 
tant to maximize a certain physical property (such as T c ) of a conventional superconductor. 
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Such a concept is very appealing and it suggests an expansion to the concept of an optimal 
spectrum which is again a delta peak spectrum with a delta peak of strength A at some po- 
sition uj*(fi*) both chosen to reproduce all known properties of a superconductor optimally. 
Such a spectrum will then provide information on the phonon mode most important for a 
specific superconductor if an a 2 F(uj) cannot be derived from experiment. It can also help 
to develop an a 2 F(u) in all cases where the phonon density of states G(u) is known. 

We would like to put this concept to test using the borocarbides LuNi 2 B 2 C and YNi 2 B 2 C 
for which extensive experimental data exist@ and for which G(uS) is known from theoretical 
work.l'i From experimental data of the upper critical field H C 2 which displays a pronounced 
upward curvature close to T c in single crystal! and polycrystallinei samples we also as- 
sume these systems to be anisotropic.0 Shulga et a/.0 explained this upward curvature of 
H C 2{T) close to T c by considering two bands, one of which being more deeply involved in 
the transport properties of the compound. The authors utilized an s-wave electron-phonon 
Eliashberg formalism and there is growing evidence that the order parameter in YNi 2 B 2 C is 
indeed of s-wave symmetry.^ It is interesting to note in passing that concept introduced by 
Prohammer and SchachingerS is effectively a two-band model described by an anisotropic 
electron-phonon interaction spectral density.00 

Sec. II of this paper reviews the theoretical background, Sec. Ill discusses the results of 
our analysis, and, finally, in Sec. IV our conclusions are drawn. 

II. THEORY 

The theoretical approach towards a theory of anisotropic polycrystalline superconductors 
within the framework of Eliashberg theory is based on the separable model for the anisotropic 
electron-phonon interaction introduced by Markovitz and Kadanoffl which was extended 
by Daams and Carbotte@ to describe an anisotropic electron-phonon interaction spectral 
function: 

a?F{u)w = (1 + a k )a 2 F(u)(l + ay), (3) 
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where k and k' are the incoming and outgoing quasi-particle momentum vectors in the 
electron-phonon scattering process and a k is an anisotropy function with the important 
feature (a k ) = 0, where (• ■ ■) denotes the Fermi surface average. As anisotropy effects are 
generally assumed to be rather small, it is sufficient to keep the mean square anisotropy (a 2 ) 
as the important anisotropy parameter. Finally, a 2 F{uj) is the electron-phonon interaction 
spectral density of the equivalent isotropic system. 

Thermodynamic properties of a superconductor are calculated from the free energy dif- 
ference AF between the normal and superconducting stated 

af = nTN (o) f; / (V*£K) + ££K) - |£ k K)|) [ 1 - -=J^^L= ] \ , (4) 
n \V )\ JcoiM + aim)/ 



with the quasiparticle density of states N(0) at the Fermi level, the renormalized quasi- 
particle frequencies u; k (c<j n ) and the Matsubara gaps A k (c<j n ) which are the solutions of the 
nonlinear s-wave Eliashberg equations: 



Wk(w„) = u n + ttT^2 ( I A k>k /(m -n) + 5 m , n -^ ) Uk ) (5a) 



A k (u; n ) =tcTY^( \Xk,k'(m - n) - /i kk , + 5^-^- J Ak \ . (5b) 
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The u; k (u; n ) are the normal state quasiparticle frequencies determined by 

^k(^n) = + ttT^ / A kjk /(m -n) + 5 m .n-^jf) sgnu; m . (6) 

In these equations u c , the cutoff frequency, is usually an integer multiple of the Debye 
frequency of the system, u n = nT(2n + l),n = 0, ±1, ±2, . . ., t kk , = l/(27r(r 4r ) kik /) is the 
anisotropic scattering rate due to inelastic impurity scattering with (rt r ) kjk ' as the anisotropic 
transport relaxation time, /x k k , is the anisotropic Coulomb pseudopotential, and 



KUm-n)=2ldn Q2+ (7) 
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In case of weak anisotropy effects the k, k' dependence of the Coulomb pseudopotential 
and of the impurity scattering is neglected and the anisotropy of the Matsubara gaps is 
described by the ansatz 



A k (u; n ) = A (u n ) + a k Ai(c^ n ), (8) 

with Ao,i(k> n ) being isotropic functions. In applying equation @ to Eqs. (|5|) only terms of 
the order of (a 2 ) are kept. 

The upper critical field H c2 (T) of an anisotropic polycrystalline superconductor employs, 
in addition, a separable ansatz to describe the anisotropy of the Fermi velocity t>F,ii^ 

v FM = (l + b k )(v F ), (9) 

with (v F ) the isotropic Fermi velocity. 5 k is an anisotropy function defined in the same way 
as a k . Again, only terms of the order (b 2 ) are kept in case of small anisotropy effects. The 
upper critical field in its temperature dependence is then described by the following set of 
equations:^ 

A k (w n ) = tcTJ2(1 + a k )A(m - n) ({I + a k /)A k /(u; TO )xk'(W)) 

m 

- nT Y, (V - i,™^) (Ak>{uj m )xk>{m))' , (10a) 



Xk(n) = —= / dxe x tan jj- — - , (10b) 







and 



« k = ^ c2 (T)(l + 6 k ) 2 K) 2 . (11) 



III. DATA ANALYSIS 

A. The a 2 F(uj) spectra 

As there are no data available which would allow to determine the electron-phonon 
interaction spectral density a 2 F(u) directly by inversion, we have to start our analysis using 
the phonon density of states G(u) which is known from theoretical work.!'! We construct 
the electron-phonon interaction spectral density using the ansatz 

a 2 F{u) = cuj s G{uj), (12) 
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with constants c and s. The coupling strength between electrons and phonons is generally 
assumed to become weaker with increasing energy u which leaves s to be negative and it is 
set to —1/2 for simplicity, a value which was also proposed by Junod et a/.0 in their analysis 
of A15-compounds. (Gonnelli et a/.0 used instead a two step weighting function which was 
determined by a fit to the high temperature part of the resistivity of a very pure single 
crystal YNi 2 B 2 C sample.) The constant c is then used to rescale the spectrum to obtain 
the experimentally observed T c (15.45 K in the case of YNi 2 B 2 C and 16 K for LuNi 2 B 2 C) 
using the linearized versions of Eqs. (|) applied to an isotropic system. Fig. [I] presents the 
low energy part of the resulting a 2 F(u) obtained for a Coulomb pseudopotential fi* = 0.13. 
For YNi 2 B 2 C A = A(0) = 1.071 (solid squares) and for LuNi 2 B 2 C A = 1.267 (solid triangles) 
which identifies both materials as medium coupling strength superconductors. 

In the analysis presented here, the actual energy dependence of the a 2 F(uj) spectrum is of 
no importance as we are going to replace the a 2 F(uj) by an Einstein spectrum with its £-peak 
of strength A at some fixed frequency w*([j,*) with A and chosen to give the measured 

T c , the appropriate value for A, and the best possible fit to experiment. Nevertheless, these 
a 2 F{u) spectra already give a pretty good idea where to place the 5-peak because a quite 
natural choice is to place the 5-peak at the center of the area under the respective a 2 F{uj) 
spectrum. 



B. YNi 2 B 2 C, clean-limit case 

We start the procedure using data for the system YNi 2 B 2 C and begin with the temper- 
ature dependence of the thermodynamic critical field 



VoH c (T) = ^2/i AF(T), (13) 
with hq the permeability constant of vacuum, and with the deviation function 

D ^f§- (1 - f2 >< < 14 ' 

where t = T/T c . The free energy difference AF(T) is determined experimentally by a 
twofold integration of the specific heat 
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T T' 

AF(T) = - J JdT'dT"^p- (15) 
o o 

measured in magnetic fields between < fioH < 9 T.iEl 

The best agreement is found for cj*(/x*) = 17meV which is close to the center of the 
area under the a 2 F{uj) spectrum (solid squares, Fig. |l|) and a A = 1.054. The numerical 
results obtained from solving Eqs. (|5]) together with (|j) are presented in Figs. and |3| 
for H C (T) and D(t) respectively. (We notice the large error bar on the D(t) data points 
which results from the method used to determine AF(T) from experiment and from the 
necessary extrapolation of the data to T — ► to find H c (0).) Experimental data are shown 
as open squares, theoretical results are presented for the isotropic case ((a 2 ) = 0, solid line), 
(a 2 ) = 0.01 (dashed line), (a 2 ) = 0.02 (dotted line), and (a 2 ) = 0.03 (dash-dotted line). It 
is obvious that results for 0.02 < (a 2 ) < 0.03 agree best with experiment. 

In the next step data for the upper critical field H C 2(T) is used to determine the anisotropy 
parameter (b 2 ) of the Fermi velocity (vp) within the above range of (a 2 ) values. Here, the 
upward curvature of if c2 (T) at T c is used to fit (b 2 ) while (vp) sets the scale according 
to Eq. (|TT|). In order to remove ambiguities, an experimental value of the average Fermi 
velocity (vp) can be derived from the plasma frequency Q p by using 

hn p = ^47re 2 (v F )N(Q)/3. (16) 

For LuNi 2 B 2 C Vl p = 4meV0 which results in (vp) ~ 0.28 x 10 6 m/s. Application of Eqs. ([TOf) 
with (b 2 ) = 0.315 and (vp) = 0.275 x 10 6 m/s results in the dashed curve in Fig. |7] for 
(a 2 ) = 0.02 while (b 2 ) = 0.305 and (v F ) = 0.285 x 10 6 m/s are found for (a 2 ) = 0.03 (solid 
line, Fig. |j) to give an optimal fit to experiment (open squares). 

The same procedure can be applied to analyze the experimental data of LuN^E^C in 
a clean limit approach using a single peak Einstein spectrum. The peak is set at the 
energy = 15meV, again close to the center of area under the a 2 F(uj) spectrum 

(solid triangles, Fig. |l|). The results of such a calculation are presented in Table | and 
compared to the YM2B2C anisotropy parameters. 
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This proves that it is indeed possible to describe the features of YNi 2 B 2 C and LuNi 2 B 2 C 
rather well using Eliashberg theory for anisotropic s-wave superconductors and a simple 
Einstein spectrum to describe the energy dependence of the electron-phonon interaction 
spectral density. Nevertheless, the agreement at low temperatures is still not perfect (see 
insets of Figs. § and ^) which demonstrate that the low temperature data are somewhat 
overestimated by our analysis. 



C. Impurity scattering 

In reality, the sample develops some residual resistivity at low temperatures which is an 
indication of some impurity content. Thus, a clean limit analysis of experimental data as it 
was presented in the previous subsection can only be a first step which allows to put some 
margins on the various anisotropy parameters. It would then be standard procedure^ to 
load the sample under investigation in a controlled way with some impurities and to measure 
the change in T c and in the residual resistivity function of impurity concentration. 

This gives another, rather reliable estimate for the anisotropy parameter (a 2 )0 and allows to 
calculate the impurity parameter t + which enters Eqs. (d) and (|10D from the Drude relation 

t+ = PrMp_ (17) 

Such data is not available for YNi 2 B 2 C and LuNi 2 B 2 C and we have to develop a different 
strategy to achieve a more realistic simulation. We make use of the fact that the critical 
temperature of an anisotropic clean limit system (t + = 0), T c o, is always greater than 
the critical temperature of a realistic system with t + > 0. One can therefore choose a 
hypothetical value for T c q and calculate how T c decreases with increasing values of t + for 
a fixed value of (a 2 ). Results of such a calculation are shown in Fig. [| (solid line for 
T c0 = 15.5 K and (a 2 ) = 0.03, dashed line for T c0 = 15.5 K and (a 2 ) = 0.035, dotted line 
for T c0 = 15.55 K and (a 2 ) = 0.035, and dash-dotted line for T c0 = 15.6 K and (a 2 ) = 0.04). 
This defines the value t + necessary in the simulation for the realistic system to obtain the 
experimental value of the critical temperature (labeled T c , YNi 2 B 2 C in Fig. |^). 
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We also have to keep in mind that adding impurities increases H&(T) even in isotropic 
systems.0H Moreover, adding impurities 'smears out' the anisotropyEI which results in 
increasing values of H C (T) and an additional increase of H C 2(T), and, furthermore, in a less 
pronounced upward curvature of H C 2(T) close to T c i3 if (a 2 ) and (b 2 ) are kept constant in 
the calculations. Thus, one will have to compensate for adding impurities by an increase 
of the anisotropy parameters. On the other hand, the anisotropy parameters found for the 
clean limit system already seem to be rather realistic (they are close to the values given 
by Manalo et a/.i found by a more elaborate analysis) and we will therefore relay to other 
means to reestablish agreement with experiment. 

Adding a second 5-peak to the a 2 F(uj) spectrum with its position and strength cho- 
sen that, again, the best possible agreement with experiment can be established by 
changing the anisotropy parameters only minimally is suggested by functional derivatives 
SH C 2(T) /5a 2 F(u)^ In particular, functional derivatives reveal that adding spectral weight 
at high energies to the a 2 F(uj) spectrum makes it less effective for H C 2(T) at low temper- 
atures while adding spectral weight at low energies has just the opposite effect. At high 
temperatures H c2 (T) is far less sensitive to changes in the spectral weight. This is of im- 
portance as our calculations overestimate H c2 (T) at low temperatures already for the clean 
limit system in the case of YM2B2C (Fig. |J). We therefore add some spectral weight at 
higher energies and a natural choice for the energy U2 at which this second 5-peak is to be 
placed is the energy of the maximum in the a 2 F(uj) spectrum, i.e. UJ2 = 21 meV for YN12B2C 
(Fig. [I]). The strength A 2 of this second peak is chosen to be a tenth of the strength A of 
the primary 5-peak. This new 25-peak spectrum has then to be rescaled to reproduce the 
experimental value of T c for the fixed value fj,* = 0.13. This procedure results in a A = 1.04 
just marginally smaller than the A of the Einstein spectrum. 

Using this new 25-spectrum H C (T) and H C 2(T) are recalculated with the parameters 
shown in Table O. The results for H C (T) are presented in Fig. || We see that the low 
temperature values of H C (T) are still slightly overestimated by our model calculations, and 
the results for if c2 (T) (Fig. |7|) reveal that if c2 (T) is now a bit underestimated at very low 
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temperatures but otherwise the agreement is almost perfect for all sets of impurity and 
anisotropy parameters. 

Again, the same procedure can be applied to the system LuNi 2 B 2 C. We use a 25-spectrum 
with the second peak placed at the maximum in the a 2 F(uj) spectrum (solid triangles, Fig. [I]) 
u>2 = lOmeV and with its strength given by A/A 2 — 7.14 because H c2 (T) was originally 
slightly underestimated. The resulting A = 1.237 and the anisotropy parameters for best 
agreement with experimental data found for this 25-spectrum are quoted in Table [TT]. 

It is important to emphasize at this point that the agreement between theoretical predic- 
tions and experiment is equally good for all sets of parameters quoted in Table ||. Thus, only 
an experiment in which the residual resistivity of the sample is measured can finally help 
in pinning down the 'real' anisotropy parameters. Our analysis only helped to substantially 
narrow the range of realistic values of the anisotropy parameters. 

IV. SUMMARY AND CONCLUSION 

We investigated two borocarbide systems, namely YM2B2C and LuN^F^C, to prove 
whether the concept of an optimal electron-phonon interaction spectrum is applicable to 
anisotropic superconductors in general and to the borocarbides in particular. Such a concept 
seems to be very helpful if only little is known about the electron-phonon interaction spectral 
function a 2 F{oj). 

This concept is indeed applicable to anisotropic clean limit systems where an Einstein 
spectrum with its peak placed near the center of the area under a model-a 2 F(a;) spectrum 
proved sufficient to obtain an excellent agreement between theoretical predictions and ex- 
periment over the whole temperature range for thermodynamics and H c2 (T). The latter 
is particularly sensitive to anisotropy effects and details in a 2 F(u>). Nevertheless, even for 
this property excellent agreement could be achieved. This proves that the model of an 
anisotropic electron-phonon interaction spectral function a 2 F(a;)k,k / can be used to explain 
the upward curvature of H c2 {T) close to T c . 

Systems with impurities require compensation of the smearing out of anisotropy by in- 
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creasing the values of the anisotropy parameters (a 2 ) and (b 2 ) if the Einstein spectrum 
concept is to be extended even to this case. This could result in rather big values for these 
parameters and in unrealistic values of (vp) necessary to reproduce H C 2(T) on an absolute 
scale. As a way out of this problem we offer to allow the optimal spectrum to be a 26- 
spectrum with the position of the main peak determined from clean limit calculations. Our 
analysis suggests various critical temperatures T c0 for the clean-limit system which is the 
'origin' of the realistic sample. Measuring the residual resistivity of the sample under inves- 
tigation will then determine t + which in turn gives the appropriate anisotropy parameters 
(see Table 0). 

Finally, the applicability of the concept of an optimal a 2 F{uj) spectrum to the systems 
YNi 2 B 2 C and LuNi 2 B 2 C proved that both are classical s-wave electron-phonon supercon- 
ductors adequately described by an Eliashberg theory of anisotropic superconductors. 
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TABLES 

TABLE I. Anisotropy parameters for YM2B2C and LUM2B2C from the clean limit calculations, 
is in meV, T c in K, and (vf) in 10 6 m/s. 



Material 




A 


T c 




(a 2 ) 


(b 2 ) 




YNi 2 B 2 C 


17.0 


1.054 


15.45 




0.02 


0.315 


0.275 




17.0 


1.054 


15.45 




0.03 


0.305 


0.285 


LuNi 2 B 2 C 


15.0 


1.174 


16.0 




0.02 


0.255 


0.288 




15.0 


1.174 


16.0 




0.03 


0.25 


0.298 


TABLE II. Anisotropy parameters for an 


anisotropic system with impurities used to simulate 


the experimental data found for YM2B2C. t + 


is given in 


meV and (vf) 


in 10 6 m/s. 




Material 


T c o 


(a 2 ) 




t+ 




(b 2 ) 


(v F ) 


YNi 2 B 2 C 


15.5 


0.03 




0.366 




0.330 


0.288 




15.5 


0.035 




0.306 




0.325 


0.293 




15.55 


0.035 




0.684 




0.345 


0.295 




15.6 


0.04 




0.961 




0.355 


0.301 


LuNi 2 B 2 C 


16.05 


0.03 




0.533 




0.275 


0.305 




16.05 


0.035 




0.455 




0.275 


0.31 




16.1 


0.035 




1.025 




0.375 


0.318 




16.1 


0.04 




0.859 




0.295 


0.32 
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FIGURES 

FIG. 1. Spectral densities a 2 F{uj) of LUN12B2C (solid triangles) and YN12B2C (solid squares) 
in the low-energy range, rescaled to obtain the measured critical temperature T c for a fixed value 
of the Coulomb pseudopotential fj,* = 0.13 from the solutions of linearized Eliashberg equations 
fro isotropic systems. 

FIG. 2. The temperature dependence of the thermodynamic critical field H C {T) in YN12B2C, 
obtained from solutions of the Eliashberg equations (0) for an Einstein spectrum with its <5-peak 
at w*(/i*) = 17meV and for various anisotropy parameters (a 2 ), namely (a 2 ) = (isotropic case, 
solid line), (a 2 ) = 0.01 (dashed line), (a 2 ) = 0.02 (dotted line), and (a 2 ) = 0.03 (dash-dotted line). 
The open squares represent experimental data. 

FIG. 3. The thermodynamic critical field deviation function D(t) as a function of the reduced 
temperature t = T/T c in YM2B2G obtained from solutions of the Eliashberg equations (||) for an 
Einstein spectrum with its <5-peak at u;*( / u*) = 17meV and for various anisotropy parameters (a 2 ), 
namely (a 2 ) = (isotropic case, solid line), (a 2 ) = 0.01 (dashed line), (a 2 ) = 0.02 (dotted line), 
and (a 2 ) = 0.03 (dash-dotted line). The open squares represent experimental data. 

FIG. 4. The temperature dependence of the upper critical field H C 2(T) in YM2B2C obtained 



from solutions of Eqs. 10 for an Einstein spectrum with its <5-peak at w*(/U*) = 17meV and for 
various anisotropy parameters, namely (a 2 ) = 0.02, (b 2 ) = 0.315, and (vp) = 0.275 x 10 6 m/s 
(dashed line), and (a 2 ) = 0.03, (b 2 ) = 0.305, and (v F ) = 0.285 x 10 6 m/s (solid line). The 
open squares without visible error bars represent experimental data obtained from specific heat 
measurements by entropy conservation. The two low temperature data points were determined 
from resistivity measurements, with error bars indicating the width of the transition.! 
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FIG. 5. The critical temperature T c of an anisotropic superconductor as a function of t + which 
is proportional to the impurity concentration for fixed values of the anisotropy parameter (a 2 ). The 
solid line is for T c q = 15.5 K and (a 2 ) = 0.03, the dashed one for T c o = 15.5 K and (a 2 ) = 0.035, 
the dotted one for T c q = 15.55 K and (a 2 ) = 0.035, and the dash-dotted one for T c q = 15.6 K and 
(a 2 ) = 0.04. The thin straight line labeles T c and indicates the experimental value of the critical 
temperature of our YM2B2C sample. 

FIG. 6. The temperature dependence of the thermodynamic critical field H C {T) in YM2B2C, 
obtained from solutions of the Eliashberg equations (|5|) for a 25-spectrum for the various anisotropy 
parameters of Table [HI The solid line is for T c q = 15.5 K and (a 2 ) = 0.03, dashed is for T c q = 15.5 K 
and (a 2 ) = 0.035, dotted is for T c q = 15.55 K and (a 2 ) = 0.035, and dash-dotted is for T c q = 15.6 K 
and (a 2 ) = 0.04. The open squares represent experimental data. 

FIG. 7. The temperature dependence of the upper critical field H C 2{T) in YM2B2C, obtained 
from solutions of equations ( |io| ) for a 2<5-spectrum for the various anisotropy parameters of Table ||. 
The solid line is for T c q = 15.5 K and (a 2 ) = 0.03, dashed is for T c q = 15.5 K and (a 2 ) = 0.035, 
dotted is for T c0 = 15.55 K and (a 2 ) = 0.035, and dash-dotted is for T c0 = 15.6 K and (a 2 ) = 0.04. 
The open squares without visible error bars represent experimental data obtained from specific heat 
measurements by entropy conservation. The two low temperature data points were determined from 
resistivity measurements, with error bars indicating the width of the transition.! 
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